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Abstract 

Methods to quantify predictability properties of atmospheric flows are pro- 
posed. The "Extended Self Similarity" (ESS) technique, recently employed in 
turbulence data analysis, is used to characterize predictability properties at short 
and long times. We apply our methods to the low-order atmospheric model of 
Lorenz (1980). We also investigate how initialization procedures that eliminate 
gravity waves from the model dynamics influence predictability properties. 

1 Introduction 

The behavior of a system such as the atmosphere cannot be predicted with prescribed 
spatial detail beyond a certain characteristic time because initially close trajectories 
diverge with time. Since the first systematic studies of error growth, low-order models 
of the atmospheric circulation, such as the Lorenz (1963,1980) models, have been used to 
study the predictability properties of chaotic systems in general and of the atmosphere in 
particular. In spite of the their lack of realism, low order models can be still considered 
useful tools to study some of the mechanisms that limit predictability. Several different 
methods have been adopted in studies on the subject. In this paper we will follow the 
guidelines of methods employed in the theory of dynamical systems (Benzi et al. 1985; 



Benzi and Carnevale 1989) in order to study the statistics of divergence of initially close 
forecasts. 

As is well-known (Oseledec 1968), for chaotic systems the asymptotic divergence 
of initially close trajectories behave exponentially with a rate given by the maximum 
Lyapunov exponent Ai (Benettin et al. 1980). Lyapunov exponents are defined as global 
averages on the attractor of a chaotic system; however, predictability is a local property 
in the phase space. Indeed, operational forecasting experience indicates that the skill 
associated with individual forecasts may vary greatly. Thus, predictability properties of 
atmospheric models exhibit large variability in both time and phase space (Palmer et al. 
1990); that is, predictability is characterized by the presence of fluctuations. 

Another important point is the existence of a period of transient growth that can 
exceed the Lyapunov exponential growth rate (Mukougawa et al. 1991; Nicolis 1992; 
Trevisan 1993; Trevisan and Legnani 1995). Error amplification on short time scales is 
controlled by rapidly growing perturbations that are not of normal-mode form. If the 
normal modes of instability of a stationary solution are not orthogonal, perturbations 
may for some time grow faster than the most unstable normal mode and the global 
average growth rate of errors, initially in a random direction, attains the asymptotic 
value given by the first Lyapunov exponent after a finite period of time. The concept 
of transient enhanced exponential growth was introduced by Farrell (1985) and Lacarra 
and Talagrand (1988). A common application of this concept used in numerical weather 
prediction consists in finding the initial perturbation that grows fastest in a given period 
of time (ECMWF 1992; Mureau et al. 1993; Molteni et al. 1996). 

In this paper we want to address the following questions: 

• how can the fluctuations in the predictability properties of atmospheric models be 
described? 

• how can the predictability properties in the transient region of super-exponential 
growth be described? 

• what is the influence of initialization procedures on the predictability fluctuations? 

Our major goal will be to give new insights on these three (related) problems. Specif- 
ically, we will use, for the first time in atmospheric science, mathematical tools which are 
well-known in the theory of dynamical systems and in other contexts (as in turbulence 
theory) and which allow the quantitative characterization of predictability fluctuations 
at long and short times. 

The outline of the paper is as follows: in sec. 2, we discuss the mathematical back- 
ground needed to describe the statistics of error growth. A quantitative definition of its 
statistical fluctuations (a characteristic also known as intermittency) is given in terms of 
an infinite set of characteristic time scales, related to the so-called generalized Lyapunov 
exponents L(q). In the absence of fluctuations, the generalized Lyapunov exponents 
are linearly dependent on Ai, thus making the maximum Lyapunov exponent the only 
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relevant parameter defining the predictability. We also discuss how to characterize the 
predictability fluctuations at short times by means of new techniques recently applied 
in turbulence theory. In sec. 3, we apply these ideas to the low-order model introduced 
and studied by Lorenz (1980). In sec. 4 we investigate how predictability properties 
are modified by the superbalance equations introduced by Lorenz (1980), which do not 
not support gravity wave oscillations. This allows us to make an assessment of the 
role played by gravity wave oscillations in enhancing predictability fluctuations at short 
times. Conclusions follow in sec. ^. 



2 Predictability fluctuations: the general theory 

In this section we shall review some basic mathematical definitions and concepts re- 
lated to sensitive dependence of the trajectory of a generic dynamical system on initial 
conditions. 

Let us consider an N-dimensional dynamical system given by the equations: 



x = f(x), (1) 

where x e R N . The time evolution of an infinitesimal disturbance 5x(i) is given by the 
linearized equations: 



5±(t) = Df 6x(t), 



(2) 



where <5x(0) ^ and Df i 



dXj 



x=x(i) 



We can define the response function R(t, 0) as 



*M)-^jf. (3) 



The Oseledec theorem (Oseledec 1968) tells us that, for t — > oo and for almost all (in 
the sense of measure theory) initial conditions 5x(0), we have: 



i?(t,0)^e Al ', (4) 

where Ai is called the maximum Lyapunov exponent (Benettin et al. 1980). One may 
venture to guess that this exponential growth is valid even at finite times and that A^ 1 
is the only characteristic time scale of the error growth; however as is well-known, this 
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is not the case (see, for example, Trevisan and Legnani 1995). Indeed, even if t is large 
enough to allow the observation of an exponential growth rate of R(t, 0), A^ 1 is not the 
only relevant time scale: depending on the initial conditions, dynamical systems can 
show states or configurations that can be predicted for times longer or shorter than A^ 1 . 

The question we want to address is related to the existence of well-defined mathemat- 
ical quantities able to characterize these predictability fluctuations in chaotic systems 
(Eckmann and Procaccia 1986; Paladin et al. 1986) and in particular in atmospheric 
models (Benzi and Carnevale 1989; Trevisan and Legnani 1995). This aspect will be 
investigated in sec. |2T1". 



Another issue to address is that for times not long enough, the dynamical behavior 
of R(t, 0) is not necessarily characterized by an exponential growth rate. This question 
has been investigated by Mukougawa et al. (1991), Nicolis (1992), and Trevisan (1993), 
Trevisan and Legnani (1995), amongst others, and turns out to be of great importance in 



realistic applications of predictability theory. We shall come to this question in sec. |2~2. 



2.1 Characterization of predictability fluctuations at long times 

In this section we recall the concept of generalized Lyapunov exponents, already intro- 
duced in atmospheric applications by Benzi and Carnevale (1989). We point out that 
Ai is an average quantity and we define the probability Pti^)) of having a local exponent 
7 different from Ai. Then we derive the relationships between all these quantities and 
we give two examples of simple probability density functions (p.d.f.'s). In particular, 
we show that the log-Poisson p.d.f, proposed here for the first time in a meteorological 
context, has some interesting features regarding the description of strong fluctuations. 

2.1.1 Basic concepts 

In this section we review some general statistical properties of the response function 
R(t, 0) (see also Benzi and Carnevale 1989). Let us consider times long enough to 
have an exponential growth rate for the response function R(t,0). Let us introduce the 
quantities: 



(R(t,0Y), (5) 

where (•) is the average over different initial conditions. 

If A i were the only relevant time scale characterizing the error growth, then we should 
expect: 



(R{t, 0) q ) ~ e Xiqt . (6) 
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On the other hand, if many times scales characterize the error growth, we should 
have the more general behavior: 



(R{t,0) q ) ~e L{q) \ (7) 

where the L(q)'s are the so-called generalized Lyapunov exponents (Fujisaka 1983; Benzi 
et al. 1985). In general, L(q) is not a linear function of q. 

Before explaining in detail the physical meaning of eq. ([7p, let us give simple examples 
of dynamical systems that satisfy either eq. @ or eq. (0) . 

Let us consider the one dimensional "tent" map (see Fig. 1): 



(0 < x < c) 

Xn+l - i ^ ( c < x < 1) 



In this case, t = n is a discrete time and the corresponding response function R(n, 0) is 
given by: 



R(n,0)=Y[Di, (9) 

i=\ 

where Di is given either by 1/c or 1/(1— c), depending on whether X{ G [0, c] or Xi G [c, 1], 
respectively. 

In order to compute (R(n, 0) 9 ), we need to know the p.d.f. of x. For the tent map 
the p.d.f. is that of the uniform distribution on the interval [0, 1] (Frisch 1995). It then 
follows from (|9|) that: 



(R(n,0y) ={j^dx + j^dxr (1Q) 
= {c l ~ q + (1 -c) l ~ q } n = e Liq)n , 1 ; 



where 



L(q) = ln[c 1 - 9 + (l-c) 1 - 9 ]. (11) 

First of all, notice that for c = 1/2, all points of the dynamical system are character- 
ized by the same slope D n = 2. In this case we have: 



L(q) = In 



1-9 



qlnl 



(12) 



Secondly, for c > 1/2 (and in particular for c very close to 1) the system is characterized 
by two different slopes, 1/c < 2 and 1/(1 — c) > 2; this means that there are states of 
the system for which the error growth is 'slow' (x G [0, c}) and states for which the error 
growth is 'fast' (x G [c, 1]). In this case, L(q) is no longer characterized by a single time 
scale. 

Coming back to the interpretation of eq. (0) and assuming that R(t, 0) already shows 
an exponential behavior, we can write: 



J R(t,0)~e 7 *, (13) 

where 7 is the local error growth exponent; the error growth exponent is local in the 
sense that it depends on the particular initial condition under consideration. We now 
briefly review how Ai can be seen as the average over all possible local exponents of the 
system. We can rewrite R(t, 0) as: 



M 

R(t,0) = Y[R(U,ti-i), (14) 

1=1 

where = t < ti < ... < t { < < ... < t M = t, and we are regarding the trajectory 
as a sequence of M trajectories. By introducing the notation: 



RiU^r) ~ e^-^\ (15) 

we have: 



e 7 * = e £<!i7i(ti-*«-i)_ (15) 
Therefore, assuming U — £j_i = At for any i, we have t = MAt and 



i=l 

Equation fllTD tells us that 7 is just the average value of ji's along the trajectory. By the 
Oseledec (1968) theorem we know that, for t — > 00 (i.e., for M — > 00), 7 — > X t . Thus, Ai 
is the average over all possible local exponents 7^ of the system. We can now ask what 
is the probability of having 7 7^ Ai for finite times. In order to answer this question, 
we introduce the probability Pti'j) of the local error growth being equal to 7 at time 
t = MAt. The introduction of a p.d.f. Pt(j) is possible under the ergodicity assumption 



6 



(for a review of the ergodic theory, see for example Halmos (1956)). The p.d.f. is then 
related to the existence of an attractor set, described by an invariant measure, which 
can be operatively defined by computing the average fraction of the time spent by the 
evolving system in any portion of the attractor. The quantity (R(t, 0) q ) can be expressed 
by the integral: 



(R(t,oy)= / P t ( 1 )e" lMAt d r (18) 



As already pointed out, only for M — > oo do we have Pt(j) — » £(7 ~~ -^l)- F° r finite 
but large M, the large deviations theory (for a general presentation see Varadhan (1984) 
and Ellis (1985); for a more physical treatment see Frisch (1995)) characterizes Ptipf) as: 



Pt{l) ~ e- s ^ MA \ (19) 

where S(j) is a concave function such that S(j) > and ^(7) = for 7 = Ai- The 
large deviations theory holds, in its simplest form, for independent random variables. 
Like the law of large numbers and the central limit theorem, the large deviations theory 
has extensions to random variables with correlations, when these correlations decrease 
sufficiently fast. These are the conditions under which Pt( 7 ) takes the form (|T9|) . In the 
language of large deviations theory, 5(7) is called Cramer function or Cramer entropy 
(Mandelbrot 1991). Different systems will have different p.d.f.'s (i.e., different 5(7)). In 
the next part of this subsection we will show how the explicit form of the L(q) exponents 
can be obtained from the Cramer function ^(7) of the system. We will also show the 
general relation between L(q) and Ai. 

Inserting ( |T9"D into (|T8|), by saddle-point integration one obtains: 

(R(t, 0)") ~ J e hi-s{i)WM dl „ e L(,)MAt ) (20) 

where 



L(q) = sup[g 7 - S( 7 )] (21) 



7 



and supjlf^)} means the smallest upper bound of the function /( 7 ). 

Since 5(7) is a concave function, the maximum value of the convex function qj — S^) 
is unique and attained at the value 7 9 at which 

~ 5(7)) , n dS( 7 ) , (99 , 



7 



Thus, from ( plj) one obtains: 



L(q) = q lq - S( lq ) (23) 

and 

^ ,24) 

aq aq aq 

Since L(0) = 0, it follows from (|23"D that £(70) = 0. Thus, as an immediate conse- 
quence of the general properties of 5(7) , one obtains 70 = Ai. From J2l|) the following 
expression for the maximum Lyapunov exponent can be obtained: 

Ai = ^| • (25) 



dq 



9 =o 



The quantities 7 g are the characteristic exponents describing the predictability fluc- 
tuations of the dynamical system: if L(q) is a linear function of q, they reduce to a single 
exponent, namely Ai (see eq. (p4|)). If L(q) does not follow the linear relationship Xiq, 
then the complete set of exponents 7 g is needed to characterize predictability fluctua- 
tions. In this case the system is intermittent. The functional form of L(q) depends on 
the system, and specific assumptions on the p.d.f. Pt{l) have to be made in order to 
perform a quantitative analysis. 

2.1.2 The log-normal case 

There is no general theory about the shape of the S(^y) function, as different systems may 
have different predictability fluctuations. If correlations are weak enough, central limit 
arguments can be applied and a Gaussian law for Pj(7) may be assumed (Paladin and 
Vulpiani 1987). Several numerical studies (see, for example, Benzi et al. 1985; Benzi and 
Carnevale 1989) show to what degree the log-normal hypothesis for the statistics of the 
error growth can reproduce the computed L(q)'s for simple dynamical and atmospheric 
models. Let us recall here the main points which follow from this hypothesis, assuming 
a quadratic behavior for 5(7): 

S(i) = (7 - Ai) 2 /2//. (26) 
From eqs. ( |19"D and (|26|) it follows that the normalized p.d.f. for 7 reads: 



p,w = vkr^- (27) 
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From the transformation rule Pt(R) = Pt(l(R)) \jr\> the log-normal distribution for 
the response function R(t, 0) follows as: 

p ' (fl) = iom e '^- (28) 

The probability distribution is thus fully characterized by two parameters only: 

X 1 = (lni^O))/*, 

/* = [((lni?(t,0)) 2 )-(lni?(t,0)) 2 ]A, 1 J 

where Ai is the maximum Lyapunov exponent and fi is the second cumulant. 

In the general case, the parameters Ai and /i, as defined in (p9|), continue to be 
important, although they may not completely represent the p.d.f.: they give respectively 
the mean value and the variance of the 7-fluctuations. For this reason, the [i parameter 
is usually referred to as the intermittency of the system. The ratio — \ may be 
taken as the borderline between weak and strong intermittency. About the latter regime, 
an important point has to be stressed. Consider the most probable response function 
value R (obtained by solving dP f^ = 0) and the mean value (R). In the log-normal case 
their expressions read: 



R = e^ 1 -^) and (#) = e Alt(1+ " /(2Al)) . (30) 

From eqs. (|30| ) it follows that a rough estimate of the response that, for example, takes 
the most probable value R as representative of the distribution, fully breaks down for 
/i/Ai > 1. This approximation predicts a decreasing error for long times (lim^oo R = 0) 
instead of the chaotic error growth characterized by a positive exponent (lim t _^ 00 (_R) = 
00). 

Using S(j) given by (p6|) in the expression (p2|) we obtain: 

7q = Ai+/ig (31) 

which, substituted into (0), gives the generalized Lyapunov exponents for the log- normal 

case: 



L(q) = X iq +^q 2 . (32) 

As a reasonable definition of the characteristic predictability time T of the system, 
one can consider the following, which takes into account the effects of fluctuations: 

T ~ ~ l - (33) 

L(l) A 1 + (l/2)/x 
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It is reasonable to assume that, in the general case, L(q) would be bounded between 
the linear shape (i.e., absence of fluctuations) and the quadratic shape (i.e., strong 
fluctuations). In many cases, the log-normal shape is a good approximation for small 
deviations of 7 around its mean value, i.e., it reproduces quite well the smallest moments 
of the distribution. On the other hand, for large q's, the moments cannot be of the log- 
normal type (Paladin and Vulpiani 1987), since log-normal moments grow more than 
exponentially with q. This implies that lim^oo 7^ is not finite. This limit is related to 
the fastest error growth R* in the system, which must be finite for all physical systems. 
The log-normal approximation fails to reproduce the tails of the distribution P t (R), that 
is, the largest fluctuations. 



2.1.3 The log-Poisson case 

The physical constraint that the fastest error growth must be finite forces one to look 
for parametrizations of Pt(j) other than Gaussian. For this purpose, let us introduce, 
for the first time to our knowledge in predictability theory, recent ideas developed in 
statistical treatment of fully developed turbulence (the basic concepts can be found in 
She and Waymire (1995)). 

One way to consider alternative probability distributions is by observing that the 
quantity: 

R{ht) ~W{tj\ (34) 

obeys the multiplicative rule: 



R(t',t) = R(t',t")R(t",t), (35) 

for any t" . Thus, one is led to consider all possible probability functions .Ft (7) which are 
left functionally invariant by the multiplicative transformation fl35l). We shall call this 
class of probability functions covariant. 

By making the (strong) assumption of weak correlation between R(t',t") and R(t",t), 
an important class of distributions turns out to be covariant: this is the class of the 
infinitively divisible distributions (IDD). The normal and the Poisson distributions are 
the most popular examples of IDD (for a comprehensive list and demonstrations of the 
properties of IDD, see for example Doob (1990)). 

In the following we shall consider the Poisson distribution as the simplest example of 
a non-Gaussian IDD that gives a more suitable description of strong fluctuations than the 
Gaussian distribution. As we shall see, it permits also to satisfy the physical constraint 
that the fastest error growth has to be finite. 
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Notice that this simple (and discrete) stochastic model is not intended to represent 
the p.d.f. of a generic (continuous) physical system. Our aim here is twofold: on one 
hand, we show that the Poisson distribution, in spite of its simplicity, shares some 
important properties with realistic p.d.f. 's; on the other hand, we show that statistical 
quantities can be successfully fitted by using the Poisson formulation. In this sense, the 
Poisson assumption gives a simple way to characterize and quantify some predictability 
properties. 

Let us rewrite the response function in terms of a Poisson random variable x, factor- 
izing its exponential time dependence: 

R(t,0) ~e at f3 x , (36) 
where (3 < 1 and x is a random variable that follows the Poisson distribution: 

P t (x = n) = ^f^. (37) 
Using eq. fl37|), the qth-order moment of the response function R can be easily calculated: 

(fiq) _ ^ ( &t )" e ^ aqtpqn _ e (aq-b)t ^ W)" _ ^ [aq _ b{1 _ pq)]t ^ 
n n 

and the expression for the generalized Lyapunov exponents becomes: 

L(q) = aq-b(l- (3 q ). (39) 
From eqs. ( ff3| ) and fl25|) we also obtain 

lq = a + b(3 q In f3 and \ 1 = a + b\nf3. (40) 



Notice that for large g's the generalized Lyapunov exponents given by fl39|) behave as 
L(q) ~ aq — b. The a parameter is associated with the fastest error growth in the system, 
R* = e at , which is finite. Thus, the log-Poisson p.d.f. does not show the pathologies 
typical of the log-normal distribution. 

A comparison between R* and (R) can give a measure of intermittency, namely: 

lni?* a 

MV = W) = ~ a - (41) 

To summarize, the considerations outlined in this section allow us to investigate 
predictability properties at large t (i.e., on the attractor set) by means of quantities such 
as Ai, /i and L(q). In particular, one can consider the ratio /i/Ai (particularly important 
if a log-normal hypothesis is made) or the ratio a (in the case of log-Poisson statistics) 
as measures of predictability fluctuations. 
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2.2 Characterization of predictability fluctuations at short times 



In this section we show, with the help of a simple example, how it is possible to extract 
quantitative information about the predictability fluctuations of a system in the super- 
exponential error growth regime. Here and for the first time in a meteorological context, 
we use recent techniques employed in turbulence data analysis, which lead to a general- 
ization of the concepts previously reported for the asymptotic regime. Again, the simple 
Poisson assumption can be made and useful statistical quantities can be evaluated. 

We are here concerned with the response function for a finite time interval. A state- 
ment about the choice of initial conditions has thus to be made. In the following, we 
make the familiar assumption of an homogeneous and isotropic probability distribution 
of the disturbance vector in phase space, which was assumed originally by Lorenz (1965). 

2.2.1 An example 

In order to discuss the predictability fluctuations at short times, we will consider the 
Lorenz (1963) system (hereafter L63). For r ~ 24.74 (the critical value for the fixed 
points to be unstable) a linear behavior of L(q) is found (see Benzi et al. 1985): the 
system does not show intermittency in the long times regime. The same behavior does 
not hold for short times, which are relevant in meteorological applications. 

This point can be detected in Fig. 2, where the second and sixth order moments of the 
response function R are shown, as an example, as a function of time for r = 28. In both 
graphs, two temporal ranges are clearly detectable (the noise in the signal is simply due 
to lack in statistics). For t > r = 1.25 there is evidence of an exponential error growth, 
whose exponent, in these log-linear plots, is given by the slope of the straight line fitting 
the data: the slopes are L(2) ~ 2Ai ~ 1.6 (see Fig. 2(a)) and L(6) ~ 6Ai = 4.8 (see 
Fig. 2(b)). These values confirm the linear law L(q) = q\i for the generalized Lyapunov 
exponents found by Benzi et al. (1985). 

At shorter times (t < r), the error grows more rapidly. Furthermore, a zoom over 
t < 0.2 reveals a non-exponential behavior of the response function. This feature makes 
it difficult to extract informations from the data inside the range < t < r. 

Notice that the separation between the two temporal regimes takes place at a time 
scale r which can be identified to be r ~ 1/Ai. This seems to be a non-trivial result, 
even if more examples based on other models are necessary to get the confidence that 
this property is general and not merely a coincidence. 

In order to determine the statistical properties of the error growth at small t, we try 
to consider different representations of our data. Let us consider, as an example, the 
log-log plot of the sixth order moment against the second order moment (Fig. 3). As 
one can see, such representation reduces the spreading of points with respect to Fig. 2, 
thus improving the quality of the linear fit. Furthermore, a zoom corresponding to the 
temporal range t < 0.2 (inside box in Fig. 3) reveals the following key point: a linear 
behavior now takes place from very early time and a measure of intermittency now 
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becomes possible. The best fit for long times (tl/Ai) gives the ratio L(6)/L(2) ~ 3 (i.e., 
no intermittency is present in the long-time regime). For t < 1/Ai this ratio gives the 
values 4.9, indicating a nonlinear behavior of L(q) and thus the intermittent nature of 
the system for short times. 

Fig. 3 represents one of the main results obtained in this paper. 

From this example we argue that it is possible to achieve a quantitative measure of 
predictability fluctuations even at short time, when non-exponential growth rates are 
observed. 

2.2.2 The ESS analysis at short times 

The results of Fig. 3 suggest to try a generalization of the concepts reported in sec. 2.1. 
Accordingly, one can continue to assume the covariance property (|35|) to be valid even 
for a time-dependence of the response R(t, 0) different from the exponential one. 

Thus, the hypothesis of IDD's governing the error growth statistics continue to hold 
in the transient regime. As an example, let us consider again the Poisson case. The 
temporal dependence e* may be replaced by a more general function g(t), such that 
g(0) = 1 and lim^oo g{t) = e t . Relation ( |36|) becomes: 



R(t,0)~g(t) a p*, (42) 
where a; is a random variable with a Poisson distribution: 



, [6 \ng(t)] n e- b ln;? W , . 

Pt(x = n)= [ - . 43 

n! 

From (^) and (0), it follows that the statistical moments of the response function 
at short t take the form: 



(R(t,0y) ~ e L(9)ln9W , (44) 

with L(q) given by (|39|) . 

The results of Fig. 3 suggest that one may extract information about predictability 
fluctuations in the system by considering ratios like \n(R g ) / ln(R p ), which allow us to 
eliminate the functional dependence on g(t). This is the basic idea of the Extended 
Self Similarity (ESS) techniques recently employed in fully developed turbulence (Benzi 
et al. 1993; Benzi et al. 1995): they concern the possibility of observing a scaling 
behavior of turbulent velocity and energy dissipation fields over an enlarged range of 
spatial scales and in systems and conditions where these scaling properties are not evident 
with standard techniques. 
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The above simple considerations allow extending the notion of generalized Lyapunov 
exponents in a simple way. Accordingly, we can define the quantities: 



ln(R) ~ L(l)lng(t)- L{q) - (45) 

The constraint L(l) = 1 leaves only two free parameters, a = a/L(l) and b = b/L{l). 
As already noted (see Eq. ([41]) ), the quantity a has a straightforward meaning and can 
be considered clS db useful intermittency indicator. 

By measuring L(q) exponents, it is thus possible to obtain quantitative information 
about predictability properties even at short times. 

Let us conclude this section by introducing another new indicator that we cinsider 
useful in studying 

predictability fluctuations at short times. We define the quantity P(t) as: 



P(t) = / P t (R)dR, (46) 



R<1 

where Pt(R) is the p.d.f. of the response function R(t, 0). 

The quantity P(t) is the probability of observing an error 5x(t) smaller or equal to 
the initial error 5x(0). By construction, we have -P(O) = 1, while in the limit of large 
t, P(t) —>■ 0. If there are no fluctuations in the predictability exponents, one obtains 
P(t) = for every t ^ 0. This is not the case for intermittent systems. 
Let us compare the cases of log-normal and log-Poisson distributions. For a log-normal 
distribution, P(t) decreases monotonically toward zero. An easy way to see it is the 
following: after substituting eq.Q2E|) into eq. ([!?]) and definining the variable 

InR — \vt , , 

z = ^=^> ( 47 ) 



one obtains the following expression of eq.(^) for the log- normal case: 

P{t) = — dz. (46 



7T 

The dependence on t is now only in the upper bound of the integral: as time increases, 
the range of integration decreases, thus making P(t) a monotonically decreasing function 
oft. 

On the other hand, the short-time behavior in the log-Poisson case is completely different: 
P(t) grows as (1 — e blng ^), up to a time t such that g(t) = /5 _1//a . Then, for large t, it 
decays toward zero, following a typical step behavior induced by the discrete character 
of the Poisson distribution. 
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The above considerations imply that in an intermittent system the probability of 
observing an error smaller or equal to the initial error can be finite. Moreover, the 
behavior of P(t) for small t is completely different in the two cases of log-Poisson and 
log-normal distributions. This effect is due to the finite maximum growth rate of the 
Poisson case: the finiteness of the parameter a allows a growth of P(t) for t <i. 

These considerations illustrate the straightforward physical meaning of the quantity 
P(t) in predictability studies concerning the short time range. Moreover, a system- 
atic, quantitative investigation may be performed at any time by evaluating the L(q) 
exponents. 



3 Application to a low-order primitive equation model 

In order to apply the theory outlined in the previous section, we have considered the low- 
order model for atmospheric flow introduced by Lorenz (1980). This model represents a 
suitable truncation of simplified primitive equations (hereafter PE) on /3-channel. As in 
Lorenz (1980), the model equations are: 



CLT 



aibiXjXk - c(di - a k )xjy k + cia-, - a^yjXk (49) 
-2c 2 y0 k - u ai 2 Xi + a^yi - z { ) 
a~- = -a k b k Xjy k - a^y^Xk + c(a k - a^y^ (50) 



-cijXj — v^a^yi 



dzj 
dr 



-b k Xj(z k - h k ) - bj{zj - hj)x k + cyj(z k - h k ) (51) 

-c(zj - hj)y k + g aiXi - K aiZi + F { 

where Xj, yi and Z{ are the coefficients of the velocity potential, streamfunction and 
height fields, respectively. Here, (i,j,k) stands for any permutation over (1,2,3). For 
our studies, we have used physical parameter values which correspond to large scale 
motion in the mid-latitude atmosphere, in agreement with previous works (Lorenz 1980; 
Gent and McWilliams 1982; Vautard and Legras 1986; Lorenz 1986; Curry et al. 1995) 
The evolution equations for an infinitesimal error have been obtained by linearization 



of eqs. ( |49H5ll) . For all simulations presented in this paper, the familiar assumption of 
homogeneous and isotropic p.d.f. of the initial disturbance vector in the phase space is 
made (Lorenz 1965). The number of realizations (i.e., different trajectories) considered 
is of the order of 10 4 . For each realization of the ensemble, the initial error has been 
randomly selected on the sphere centered at the initial point of the trajectory. The 
numerical integrations have been performed using a fourth-order Runge-Kutta scheme. 

In order to show the capability of the parameter /i/Ai to capture the fluctuations in 
the predictability, let us consider the three different regimes of the PE model discussed by 
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Krishnamurthy (1985), corresponding to three different values of the forcing parameter: 
F\ = 0.1, Fi = 0.25 and F\ = 0.3. The physical regimes associated with such values 
are the following (for details, see Krishnamurthy (1985)): for Fx = 0.1 the attractor of 
the PE system is free from gravity waves; for F\ = 0.25 the trajectories vary chaotically, 
with an almost periodic behavior interrupted irregularly by short periods characterized 
by high frequency gravity waves; finally, for F\ = 0.3, gravity waves are present all the 
time. 

To evaluate the two parameters Ai and fi, defined in ( p9|) , each realization (trajectory) 
has been made to last for a time t of the order of a hundred days, large enough to apply 



all the considerations of sec. 2.1 



In Tab. 1, the values of the positive maximum Lyapunov exponent, the intermittency, 
the ratio fi/Xi, and the predictability time T for the PE system in the three aforemen- 
tioned regimes are shown. As already noted, the value of ///Ax = 1 could be considered 
as the borderline between weak and strong intermittency. We can see from Tab. 1 that 
three different values of fi/X\ occur in the three regimes. The signature of the gravity 
wave activity, in the cases F\ = 0.25 and F\ = 0.3, is emphasized by the enhancement of 
both Ai and fi, with respect to the case F\ = 0.1. The predictability of the system is thus 
reduced by persistent gravity waves. It is also worth noticing that strong intermittency 
(i.e., fi/Xi > 1) is present only in the intermediate case F\ = 0.25, characterized by the 
appearance of bursts of chaoticity containing high frequency gravity waves. 

For the regime F\ = 0.1, we have also performed a systematic numerical investigation 
of the generalized Lyapunov exponents at long and short times. The logarithm of the 
sixth order moment of the response function R as a function of time is reported in Fig. 4. 
Two regions are clearly detectable. For t large enough (t > 1/Ai ~ 16 days) the points 
can be fitted by a straight line, indicating an exponential error growth: the slope of 
the fitting line gives the sixth order generalized Lyapunov exponent L(6) (the noise in 
the signal is simply due to lack in statistics). At shorter times (t < 1/Ai), the error 
grows more rapidly; moreover, a zoom over the first 5 days reveals a strongly nonlinear 
behavior. As we shall see in the next section, this is mainly associated with the presence 
of transient gravity waves, whose amplitude is significant, especially in very the first 
times of integration. 

Notice that, as in the L63 model, the characteristic time scale associated with a 
change in the error growth regime is of the order of 1/X±. 

The use of ESS techniques allows us to improve the quality of the fit in the long- 
time region and to extract information on generalized Lyapunov exponents even at small 
t. As an example, the sixth order moment against the first order moment is shown in 
Fig. 5 in a log- log plot. We can see from the zoom over the 5 day temporal range 
that a linear behavior occurs from very early times. We have then evaluated for 14 
moments (q=0.2, 0.4, 0.6, 0.8, 1, 2,..., 10) the L(q) exponents at short and long times, 
by performing linear fits in the two temporal regions < t < 10 days and 40 < t < 300 
days, respectively. The results are shown in Fig. 6 : the points are well fitted by the 
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log-Poisson formula, with parameters 5; = 1.2, hi = 0.3 (for the long-time L(q) values) 
and a s = 4.5, b s = 4.3 (for the short-time L(q) values). In particular, notice that a s > a h 
as the comparison of the two slopes at high g's shows. This is a quantitative indication 
of the more intermittent nature of the system for short times (see eq. (pHf)). 

Focusing our attention on the short-time regime, we have also evaluated the function 
P(t) defined in sec. [2.2| , by counting the number of realizations (trajectories) having 
R(t, 0) < 1 at a certain time t: 



P(t) = ^hL. (52) 
Ntot 

Fig. 7 shows the P(t) for the PE model. Due to the appearance of transient gravity 
waves, a very strong loss of predictability occurs immediately after the beginning of the 
PE model integration: P(t) is close to and starts to grow until i ~ 5 days, a behavior 
in qualitative agreement with a log-Poisson character of the R(t,0) p.d.f.. After about 
10 days, gravity waves become lost and only slow oscillations remain. 

From the above results, it appears that the presence of gravity waves (on the at- 
tractor set as well as in the initial transient) has strong effects on predictability. In 
particular, the loss of predictability immediately after the first times of the model inte- 
gration may be relevant for practical purposes. The final question we want to address is 
thus the following: may the elimination of gravity waves produce an enhanced short-time 
predictability? We shall try to answer this question in the next section. 



4 Predictability fluctuations and their relations with 
gravity wave activity 

Our first aim is to study how the predictability fluctuations are influenced by the disap- 
pearance of gravity waves in models not supporting such waves. 

We shall focus our attention on Lorenz's (1980) algorithm exploiting the complete 
separation between quasi-geostrophic and gravity wave frequencies, and leading to the 
well-known superbalance equation, defining the slow manifold (Leith 1980; Lorenz 1980). 

Defining X G = (Xf , X 2 G , X 3 G , X 4 G , X G , X G ), where Xf = x { and X G 3 = Zi - Vi for 
i = 1, 2, 3, while X fi = (yi,y2, 1/3), Lorenz (1980) exploited the scale separation imposing 

the condition dt p k = which can be solved iteratively for each p by the Newton method 
(Lorenz 1980). 

The low-order dynamical system obtained from the PE model by calculating the 
X G components through Lorenz's algorithm (with p = 1) is called here SuperBalance 
Equation (SBE) model. 



We apply the notions outlined in sec. 2.1 and |2.2| to the PE and SBE models in order 



to quantify the effects of gravity wave activity on the predictability fluctuations. 
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The positive maximum Lyapunov exponent, the intermittency, the ratio and 
predictability time T (see eq. (|53D) are shown in Tab. 2 for the two models. Actually as 
we can see from this table, these asymptotic quantities do not permit us to discriminate 
between the PE and SBE models. The very slight differences in these chaotic indicators 
confirm the capability of SBE to reproduce the asymptotic statistical properties of the PE 
system. Let us then come back to the original question: in what sense does elimination 
of transient gravity waves make the SBE model less chaotic than the PE model? 

An answer to this question can be given by exploiting the short-time analysis pre- 
sented in sec. |2.2|. In Fig. 8 we report the same plot of Fig. 4, but for the SBE model. 



As in the PE model, two regions are clearly detectable, with the error growth for short 
times faster than that for long times. However, the zoom over the first 5 days shows an 
exponential (that is linear in the graph) error growth since the very first time, unlike the 
PE case. This is due to the projection onto the slow manifold which causes the gravity 
waves to die out completely. 

In order to improve the quality of the fit and make a comparison possible, we perform 
for the moments of the error growth the same ESS analysis already presented for the PE 
case. In particular, L(q) exponents of the SBE model in the long-time regime turn out 
to be indistinguishable from those of the PE model shown in Fig. 6 (white circles) . This 
fact confirms the success of the SBE model in reproducing the asymptotic statistical 
properties of the Lorenz model. 

Let us focus our attention on the short-time behavior. The L(q) exponents for the 
SBE model, obtained from a linear fit at short time, are reported in Fig. 9. In order to 
simplify the comparison, the corresponding curve for the PE model (already reported 
in Fig. 6 (black circles)) is shown. Several comments are worthwhile. The log-Poisson 
fit for the SBE model gives a s = 2.9 and b s = 2.4: in particular, we stress that the 
slopes a s are smaller than in the PE case. Hence, disappearance of gravity waves in the 
SBE model decreases temporal intermittency, damping predictability fluctuations in a 
sensible way. 

The enhanced predictability during the first days becomes more evident if one con- 
siders Pit) behavior for the PE and SBE models, shown in Fig. 10. As we can see, the 
SBE curve strongly differs from that of the PE system during the first days. Notice that 
the P(t) starts near 1 for the SBE model. It then follows that the very strong loss of 
predictability (i.e., P(t) close to 0) immediately after the first times of the PE model 
integration is certainly due to the transient gravity waves: for the SBE model, where 
gravity waves are not present, there is only a slight loss of predictability (i.e., P(t) close 
to 1) during the first few days. 

After about 10 days, gravity waves disappear in the PE system and the behavior of 
P(t) for the two models becomes quite similar. 

From the behavior of both L(q) and P(t), we can conclude that the SBE model system- 
atically produces at short times a more predictable flow with respect to the PE model. 
To our knowledge, this is the first statistical assessment concerning the effects introduced 
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by initialization procedures on the predictability properties of an atmospheric model. 

5 Conclusions 

In the Introduction we have posed three different questions which we believe are relevant 
for the predictability problem in atmospheric flows. 

The answer to the first question we posed (how to describe predictability fluctuations) 
can be given in terms of generalized Lyapunov exponents, L(g)'s, which represent the 
growth rate of the moments of the response to an initial error. The statistics of R(t, 0) 
are characterized by an infinite set of exponents related to (R(t,0) q ). The non linear- 
ity of L(q) vs q reflects the intermittent nature of the temporal evolution of the error 
growth. In this sense, this infinite set of exponents contains a level of information which 
is not present in the maximum Lyapunov exponent Ai, Ai being able to quantify only the 
average predictability properties of a system. The use of I/(?)'s is not completely new. 
They have been already introduced in the framework of atmospheric models by Benzi 
and Carnevale (1989). Regarding specific hypothesis on the R(t,0) p.d.f., they consid- 
ered the log-normal case. Here we have proposed the use of a log-Poisson p.d.f. (recently 
introduced in the statistical theory of turbulence for the treatment of spatial intermit- 
tency, see for example She and Waymire (1995)) to overcome some of the deficiencies of 
the log-normal p.d.f.. 

We have also focused our attention on the short-time predictability properties which 
is a rather complex problem. The fact that the error growth need not be exponen- 
tial for short times makes the analysis more difficult, introducing a complex functional 
dependence of the response function on the time t, which can vary depending on the 
system. 

Using recent ideas introduced and applied in the context of fully developed turbu- 
lence, namely the Extended Self Similarity techniques, we have shown that the concepts 
of the generalized Lyapunov exponents can be extended to the short-time regime, intro- 
ducing the L(q) exponents. They are well-defined mathematical quantities that allow us 
to quantify the intermittency in the system also at short times. This is one of the major 
results achieved in this paper. 

As a further measure of predictability fluctuations at short time, we have also con- 
sidered the probability to have a decrease of the initial error at time t (what we called 
P(t)): this is a clear indicator concerning the possibility to have convergence of two 
trajectories at finite time, due to large fluctuations around the average error growth. 

We have applied these concepts to the case of a simple atmospheric model (Lorenz, 
1980), which captures some features of synoptic scale motions at mid-latitudes. We 
chose this model because its numerical study does not need the computational resources 
necessary in the case of more complex and realistic models. Moreover, in spite of its 
simplicity, the model has a rich dynamics, characterized by appearance of the so-called 
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gravity wave modes. This property has allowed us to answer the third question posed 
in the Introduction, that is, which is the connection, if any, between short-time pre- 
dictability fluctuations and initialization procedures. To this end, we have considered 
the initialization procedure proposed by Lorenz (1980) leading to the SBE model. 
Our results concerning the PE and SBE models can be summarized as follows: 

- the generalized Lyapunov exponents evaluated for the PE and SBE models are 
well-reproduced by a statistical law of the log-Poisson type and differ both from 
the non- intermittent linear law L(q) oc q and the quadratic form implied by a 
log-normal assumption; 

- we have observed in the models two different regimes, associated with two temporal 
ranges: in the short-time regime, the error grows more rapidly and intermittently 
than in the asymptotic region. This means that standard (asymptotic) chaotic 
indicators as Ai and /x can be meaningless when considering short-time regimes; 

- comparing the asymptotic behaviors of the models, we conclude that the SBE 
procedure is successful in reproducing the long-time statistics of the original PE 
model, since all the asymptotic chaotic indicators are practically indistinguishable 
in the two cases; 

- elimination of gravity waves in the SBE model strongly influences predictability for 
short times, giving a more predictable and less intermittent flow. The initialization 
procedure is then successful in decreasing the intermittency in the short time region, 
which is strongly influenced by transient fast oscillations in the PE model; 

- strongly non-exponential error growth was observed in the case of the model sup- 
porting gravity waves (PE) after the first times of the model integrations. This 
is the signature of the strong gravity wave activity. Non-exponential behavior is 
strongly reduced in the SBE model, which does not support gravity waves. 

We would like to conclude with some remarks and suggestions for future work. It is 
necessary to better understand the physical meaning both of the transition observed in 
the predictability properties (short-time and long-time regimes) and the characteristic 
time at which this transition occurs. This characteristic time seems to be strongly related 
to the maximum Lyapunov exponent. We think that this feature (which we are currently 
investigating) is not peculia only to the models studied in this paper. 

Another point regards the study of the growth of non-infinitesimal errors: indeed, in 
this paper we have considered the error growth in the so-called tangent space, that is, 
we have studied the behavior of an infinitesimal initial error, governed by the linearized 
version of the dynamical equations. It would be interesting to perform a similar analysis 
for finite initial errors, in order to obtain more realistic information on the predictability 
when the error in the initial conditions cannot be considered infinitesimal. 
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Finally, we conclude with two remarks concerning the feasibility of our analysis when 
more complex models involving a large number of variables are considered. First, in our 
analysis we computed the moments of the response function up to the tenth order just 
to highlight the effect of predictability fluctuations. Actually, these can be captured also 
restricting the analysis to lower order moments, thus reducing the related computational 
effort. Second, since the predictability is characterized by fluctuations, the need for a 
good statistical ensemble of simulation runs is always crucial. This is independent of the 
method used to perform the statistical analysis and, in particular, is not peculiar to our 
method. The reason for this is that for systems characterized by strong intermittency, 
large fluctuations may occur with a non-vanishing probability and their effects on the 
predictability have to be properly accounted for in a statistical way. A large amount of 
statistics has thus to be gathered just in order to sample the tails of the p.d.f., where 
the rare events are placed. From this point of view, it is clear that over-simplified 
atmospheric models suggest a possible strategy to tackle the problem of predictability 
fluctuations. 



21 



Acknowledgements 

We are grateful to L. Biferale, R. Festa, C.F. Ratto, 0. Reale, M. Vergassola and A. 
Vulpiani for illuminating discussions, and Dr. Claudio Paniconi for reviewing the text. 
We thank the "Meteo-Hydrological Center of Liguria Region" where part of the numerical 
analysis was done. The author (2) acknowledge support from the Sardinian Regional 
Authorities. 



22 



References 



Benettin, G., L. Galgani, A. Giorgilli and J. M. Strelcyn, 1980: Lyapunov charac- 
teristic exponents for smooth dynamical systems; A method for computing all of them. 
Part I: Theory, and Part II: Numerical Application. Meccanica 15, 9-20 and 21-30. 

Benzi, R., G. Paladin, G. Parisi and A. Vulpiani, 1985: Characterisation of intermit- 
tency in chaotic systems. J. Phys. A, 18, 2157-2165. 

Benzi, R., and G. Carnevale, 1989: A possible Measure of Local Predictability. J. 
Atmos. Sci., 46, 3595-3598. 

Benzi, R., S. Ciliberto, R. Tripiccione, C. Baudet, C. Massaioli, and S. Succi, 1993: 
Extended self-similarity in turbulent flows. Phys. Rev. E, 48, R29-R32. 

Benzi, R., S. Ciliberto, C. Baudet, and G. R. Chavarria, 1995: On the scaling of 
three-dimensionale homogeneous and isotropic turbulence. Physica D, 80, 385-398. 

Curry, J.H., S. E. Haupt, and M. N. Limber, 1995: Low-order models, initialization 
and slow manifold. Tellus, 47A, 145-161. 

Doob, J.L., 1990: Stochastic Processes. J. Wiley and Sons, 654 pp. 

Eckmann, J. P., and I. Procaccia, 1986: Fluctuations of dynamical scaling indices in 
nonlinear systems. Phys. Rev. A, 34, 659-661. 

ECMWF, 1992: New developments in predictability. Workshop Proceedings, 13-15 
November 1991. 

Ellis, R.S., 1985: Entropy, Large Deviations and Statistical Mechanics. Springer, 364 

pp. 

Farrell, B. 1985: Transient growth of damped baroclinic waves. J. Atmos. Sci., 
42,2718-2727. 

Frisch, U., 1995: Turbulence, Cambridge Univ. Press, 296 pp. 

Fujisaka, H., 1983: Statistical dynamics generated by fluctuations of local Lyapunov 
exponents. Progr. Theor. Phys., 70, 1264-1275. 

Gent, PR., and J. C. McWilliams, 1982: Intermediate Model Solutions to the Lorenz 



23 



Equations: Strange Attractors and Other Phenomena. J. Atmos. Sci., 39, 3-13. 

Halmos, P.R., 1956: Lectures on Ergodic Theory, Chelsea, 302 pp. 

Krishnamurthy, V., 1985: The slow-manifold and the persisting gravity waves. PHD 
Thesis, MIT, 146 pp. 

Lacarra, J. and O. Talagrand, 1988: Short range evolution of small perturbations in 
a barotropic model. Tellus, 40A, 81-95. 

Leith, C.E., 1980: Nonlinear Normal Mode Initialization and Quasi-Geostrophic 
Theory. J. Atmos. Sci., 37, 958-968. 

Lorenz, E.N., 1963: Deterministic non-periodic flow. J. Atmos. Sci., 20, 130-141. 

Lorenz, E.N., 1965: A study of the predictability of a 28- variable atmospheric model. 
Tellus 17, 321-333. 

Lorenz, E.N., 1980: Attractor sets and Quasi-Geostrophic Equilibrium. J. Atmos. 
Sci., 37, 1685-1699. 

Lorenz, E.N., 1986: On the existence of a slow manifold. J. Atmos. Sci., 43, 1547- 
1557. 

Mandelbrot, B., 1991: Random multifractals: negative dimensions and the resulting 
limitations of the thermodynamic formalism. Proc. R. Soc. Lond., A 434, 79-88. 

Molteni, F., R. Buizza, T.N. Palmer, and T. Petroliagis, 1996: The ECMWF Ensem- 
ble Prediction System: Methodology and validation. Q.J.R. Meteor. Soc, 122, 73-119. 

Mukougawa, H., M. Kimoto, and S. Yoden, 1991: A relationship between local error 
growth and quasi-geostrophic states: case study in the Lorenz system. J. Atmos. Sci., 
48, 1231-1237. 

Mureau, R., F. Molteni, and T.N. Palmer, 1993: Ensemble prediction using dynam- 
ically conditioned perturbations. Q.J.R. Meteor. Soc, 119, 299-323. 

Nicolis, C, 1992: Probabilistic aspects of error growth in atmospheric dynamics. 
Quart. J. Roy. Meteor. Soc, 118, 553-568. 

Oseledec, V.I., 1968: A multiplicative ergodic theorem. Lyapunov characteristic 



24 



numbers for dynamical systems. Trans. Moscow Math. Soc, 19, 197-231. 



Paladin, G., L. Peliti and A. Vulpiani, 1986: Intermittency as multifractality in his- 
tory space. J. Phys. A, 19, L991-L996. 

Paladin, G., and A. Vulpiani, 1987: Anomalous scaling laws in multifractal object. 
Phys. Rep., 156, 147-225. 

Palmer, T.N., R. Mureau and F. Molteni, 1990: The Montecarlo forecast. Weather, 
45, 198-207. 

She, Z.S., and E. C. Waymire, 1995: Quantized energy cascade and log-Poisson 
statistics in fully developed turbulence. Phys. Rev. Lett., 74, 262-265. 

Trevisan A., 1993: Impact of transient error growth on global average predictability 
measures. J. Atmos. Sci., 50, 1016-1028. 

Trevisan A., and R. Legnani, 1995: Transient error growth and local predictability: 
a study in the Lorenz system. Tellus, 47A, 103-117. 

Varadhan, S.R.S., 1984: An Introduction to the Theory of Large Deviations. Springer- 
Verlag, 196 pp. 

Vautard, R., and B. Legras, 1986: Invariant Manifold, Quasi-Geostrophy and Ini- 
tialization. J. Atmos. Sci., 43, 565-584. 



25 



c 

X 



Figure 1: The one-dimensional "tent map". 
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Figure 2: The logarithms of the second (a) and sixth (b) order moments of the response 
function R(t, 0) as functions of time t for the L63 model. Dotted lines result from linear 
fits of the data for t > \/X\. their slopes are ~ 2Ai and ~ 6Ai, respectively. Inside boxes 
are zooms over the very short-time behaviors. 
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Figure 3: Log-log plot of the sixth order moment against the second order moment of the 
response function R(t, 0) for the L63 model. Dashed lines result from linear fits of the 
data in the short and long time regimes. Inside box is a zoom over the very short-time 
behavior. 
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Figure 4: The logarithm of the sixth order moment of the response function R(t, 0) as a 
function of time t for the PE model. Inside box is a zoom on the first 5 days behavior. 
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Figure 5: Log-log plot of the sixth order moment against the first order moment of the 
response function R(t, 0) for the PE model. Dashed lines result from linear fits of the 
data in the short and long time regimes. Inside box is a zoom on the first 5 days behavior. 
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Figure 6: The L(q) exponents as a function of q for the PE model: black circles are the 
values for < t < 10 days, while white circles denote the long-time values (40 < t < 
300 days). 
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Figure 7: Time behavior of P(t) for the PE model. 
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Figure 8: The same as in Fig. 4, but for the SBE model. 
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Figure 9: The L(q) exponents as a function of q for the PE and SBE models in the 
short-time regime (0 < t < 10 days): circles and stars are the PE and SBE values, 
respectively. 
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Figure 10: Time behavior of P(t) for PE (thin curve) and SBE (dotted curve) models. 





Xi(days x ) 


[/(days l ) 


/i/Ai 


Predictability (days) 


0.1 


6.3 x 10~ 2 


2.5 x 10~ 2 


0.4 


13.2 


0.25 


1.7 x lO" 1 


2.8 x lO" 1 


1.6 


3.2 


0.3 


4.2 x 10- 1 


3.1 x 10- 1 


0.7 


1.7 



Table 1: Maximum Lyapunov exponent Ai, intermittency [/,, their ratio /i/Ai and pre- 
dictability for the PE model for different values of the forcing parameter. 



Models 


Xi(days x ) 


[/.(days r ) 


///Ai 


Predictability (days) 


PE 


6.3 x 10~ 2 


2.5 x 10~ 2 


0.4 


13.2 


SBE 


6.5 x 10~ 2 


2.0 x 10~ 2 


0.3 


13.3 



Table 2: Maximum Lyapunov exponent Ai, intermittency [/,, their ratio f//\i and pre- 
dictability for the PE and SBE models at F 1 = 0.1. 
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